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Abstract 

We study the behavior of the tree-level S-matrix of a variety of theories as two particles become 
soft. By analogy with the recently found subleading soft theorems for gravitons and gluons, we 
explore subleading terms in double soft emissions. We first consider double soft scalar emissions 
and find subleading terms that are controlled by the angular momentum operator acting on hard 
particles. The order of the subleading theorems depends on the presence or not of color structures. 
Next we obtain a compact formula for the leading term in a double soft photon emission. The 
theories studied are a special Galileon, DBI, Einstein-MaxwelDScalar, NLSM and Yang-Mills- 
Scalar. We use the recently found CHY representation of these theories in order to give a simple 
proof of the leading order part of all these theorems. 
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I. INTRODUCTION AND SUMMARY OF RESULTS 


In 2014 subleading soft factors were found for the single soft emission of gravitons in 
Einstein gravity [1]. The structure of a (n + l)-graviton amplitude in the limit when kn+i is 
soft is given by 

•W»*1 = (41., + 4h.y + 41.,)-W» + (1) 

where r is a parameter that controls the soft limit, and 

n - 1,1- 

40) 


.5: 


n , 


( 2 ) 


gravity Z-/ jL u 

a=l ^n+1 

is Weinberg’s universal soft factor [2], In this formula is the polarization tensor of the 
soft particle. The next two soft factors were found by Strominger and one of the authors in 
[1] and are given by 


c,(l) ^ e^yk^{kn+l,pJa'') ^(2) 1 ^ (^iiu(kn+l,pJ^^)(kn+l,aJa'") 

•^gravity “ Zj h . . h ’ ‘^gravity “ o Zj h , . h 

a=l ^n+1 ^ a=l ^n+1 


( 3 ) 


In these formulas Ja^ is proportional to the total angular momentum operator acting on the 
particle. For example, when a is a scalar particle 


- k^ ^ k'' — 

a,scalar a q, a 


( 4 ) 


a,/i 


By analogy with the gravity construction, Casali identihed a sub-leading factor for the 
emission of a soft gluon in Yang-Mills [3]. For a U{N) color-ordered partial amplitude the 
structure becomes 


A4(l, 2,..., n, n + 1) = (4^4 + 2,..., n) + 0(r), 


( 1 ) 


where 


Y 


( 0 ) _ 


e ■ kr, 


e ■ k-\ 


YM 




( 1 ) _ ^p(,kn+l,pJn^') (^p{kn+l,pJi^) 


( 5 ) 


( 6 ) 


kn+1 ■ kn ^n+1 ’ k\ ^n+1 ’ ^n+1 ’ k\ 

Sub-leading soft theorems have a long history dating back to the 50’s (see e.g., [4-6]). 

Particularly known examples are the subleading terms in soft photon emissions which are 
referred to as the Low-Burneet-Kroll theorem [4]. 

Although the original proofs of (1) and (5) were performed in four dimensions, the validity 
in any number of dimensions was soon established [7-9]. In this work we focus on amplitudes 
at tree-level so the only comment on loop amplitudes is that single soft subleading terms 
have also been studied at higher orders in perturbation theory in, e.g., [10]. 
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In this paper we study two classes of theories that contain scalar particles. The hrst class 
consists of theories with neither color nor flavor structure while the second contains theories 
with a U{N) color (or flavor) group for the scalars. For both classes we consider single and 
double soft scalar limits. For the single soft limit we are only concerned with order of the 
leading term while for the double soft limit we study the actual structure of the leading and 
subleading soft factors. The leading behavior of soft scalars has been studied for a long time 
with perhaps the most well-known result being the Adler’s zero [11], and very recently were 
studied as a classihcation tool for scalar theories in [12]. Double soft limits of scalars have 
also been studied in AA = 8 supergravity as a way to explore its moduli space of vacua and 
the F^ 7 ( 7 ) structure of the theory [13]. 

In the hrst class of theories we have a special Galileon theory [14] which we will refer 
to as sGal [12, 15], the Dirac-Born-Infeld (DBI) theory, and an Einstein-Maxwell-Scalar 
(EMS) theory. Here we consider theories with a single scalar held, thus for example EMS is 
the dimensional reduction of Einstein’s theory from (D+l) to D dimensions. 

For theories in the hrst class we propose that as two scalars, say particles n + 1 and n + 2, 
become soft as ^iny (n + 2)-particle amplitude behaves as 


Mn,2 = (fcn+1 • kn,2r{S^°^ + + 0(r“) (7) 


where m = 1,0,-1 for sGal, DBI and EMS respectively, and 

(/C(j • (^kn+i ~ ^n+ 2 )) 


1 ” 

5(0) = i y 

A ^ 


1 n 


4 \ ■ (^n+1 ^n+ 2 ) ^n+1 ’ kn+2 

ka ' (^n+1 ~ ^n+ 2 ) 


1 ^ 


2 a=l ■ (,kn+l + kn+ 2 ') ^n+1 ■ kn+2 

1 


+ ka ■ {kn+\ + ^n+2) ^n+1 ■ kji+2 j 1 


2 a-1 ' (^n+1 ^n+ 2 ) ^n+1 ' ^n+2 




( 8 ) 

(9) 

( 10 ) 


Here F'G) is a multiplicative operator which has an expansion in r starting at FB) 

is a hrst-order diherential operator starting at and S'G) is a second order diherential 

operator starting at O(r^). 

Two comments are in order at this point. The hrst is that in (7) the kinematic invariant 
(fc„+i • kn+ 2 ) plays the role of a natural “dimensionful parameter” needed to link amplitudes 
with diherent number of particles. This simple dimensional argument leads to universal 
formulas for S'G) and S'F), i.e., they are theory independent. The second is that the only 
dependence on the theory under consideration appears in the multiplicative piece of FG). 
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As mentioned above, the second class of theories corresponds to those with a U{N) color 
(or flavor) structure. In this class we have the non-linear sigma model (NLSM) [16] and 
Yang-Mills-Scalar (YMS), which is the dimensional reduction of Yang-Mills theory. The 
double soft scalar emission for a color-ordered partial amplitude is proposed to be 

7W(1,2,..., n, n + 1, n + 2) = 2,..., n) + C>(r“) 


with m = 0, -1 for NLSM and YMS respectively, and 


^(0) ^ 


5(1) = 


1 ( ■ (^n+1 ^n+ 2 ) ^n+l ' 


^ kn ■ (^n+1 + ^n+ 2 ) ^n+1 ’ kji+2 

kn+l,fikn+2,u 


n+2 ^ kl ■ {kn+2 


kn+l) ^n+2 ■ k. 


n+1 


kn * (fcn+1 kn+2') kn+l ■ kn+2 


jr + 


kl ■ {kn+2 + kn+l) + kn+2 ■ kn+l 
kn+2,fj,kn+l,u 

kl • {kn+2 + kn+l) kn+2 ■ ^n+1 


jr. 


( 12 ) 

(13) 


In this formula S'(°) starts at 0{t^) while iS'(i) starts at order 0{t). Expanding 5(°) in (12) 
around r = 0 the leading order becomes 


S(o) 


1 ( kn ■ {kn+l kn+2) ^ kl • {kn+2 kn+l) j ^ C^{'f) 

2 y kn ■ {kn+l kn+2) kl • {kn+2 kn+l ) / 


(14) 


This is the famous double soft factor which although hnite depends on the relative directions 
of the soft particles. This structure is the one carrying the information of the non-linearly 
realized symmetries of the theory which is hinted by Adler’s zero in the NLSM. Although 
the form has been known for many years and tested in the known amplitudes, a proof to 
all multiplicities in the NLSM was only recently found using BCFW techniques [17]. For 
conventions and details of all the theories we consider here, please refer to [15]. 

Last but not least, we hnd the double soft photon emission theorem in DBI and EMS 
theory, with a universal leading-order soft factor. Here we record the result: when particles 
n+1 and n+2 are soft photons, an (n+2)-point amplitude becomes 


A^n+2 = {kn+l ■ kn+2)^-^S^^^Mn + 0 {t^^^^) , 


^(0) {ka ■ {kn+l kn+ 2 ))^ 


a=l 


ka ' (^n+1 ^n+ 2 ) 


Pf5,. (15) 


where m = 0, -1 for DBI and EMS respectively, and the 4x4 anti-symmetric matrix Sa is 


5a := 


^ 0 kn+l ' kn+2 k^ • €.n+l ~kn+l ' ^n+2 ^ 

~kn+l ' kn+2 0 kn+2 ' ^n+1 k^ ’ 67^+2 

’ ^n+1 ~kn+2 ' ^n+1 0 ^n+1 ' ^n+2 

\ kn+l ' ^n+2 ~ka ' ^n+2 ~^n+l ' ^n+2 0 / 


( 16 ) 
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where ik^Y := ) and (A:^')^ := A ka - fc^.i) 

are two vectors satisfying k^ ■ kn-^i = 0 and k^ ■ kn +2 = 0. Gange invariance of the soft photon 
factor becomes manifest due to the Pfaffian structure. 

In the rest of the paper we provide various levels of evidence for all the soft theorems 
proposed here. In Section IIA we list a variety of explicit checks performed with known 
formulas for the amplitudes. In Section IIC we use the recently found Cachazo-He-Yuan 
(CHY) representation [18] of the theories under study [15, 19] to give a very simple proof of 
the leading order terms in AG) for all the theories mentioned above. A proof for the sub¬ 
leading order including the leading part of A^) is also done using the CHY representation 
and is presented as supplementary material [20], where we also include the proof for the 
double soft photon theorem. We end with discussion on future directions in Section III. 


II. CHECKS AND PROOFS FOR THE NEW THEOREMS 

Here we provide strong evidence for the double soft theorems. We perform non-trivial 
checks using explicit amplitudes, and use the CHY formula to prove some of the theorems. 

As in the case of soft gluon and graviton emission, amplitudes that enter into our soft 
theorems are distributions containing momentum-conserving delta functions. For both the 
checks and proofs, it is useful to write versions of the soft relations for stripped amplitudes. 
Let us define M. = for n-point and (n+2)-point amplitudes, and recall that the 

two soft momenta are k ^^2 - 

As distributional relations, (7) and (11) are equivalent to relations for all independent 
distributions when an expansion in r is performed, i.e., the momentum-conserving delta 
functions and its derivatives. One can show that the relations for the derivatives are guar¬ 
anteed by the relations for the delta functions, or the soft theorems for stripped amplitudes 

[Mn, 2 ](r) = [(Yp • g)”^(AG) + aG) + aG)) mJ (r) + , (17) 

[M(l,2,... ,n+l,n+2)](r) = [{Yp ■ + A^)) M(l,2,... ,n)] (r) + C>(r“) ^ 

In both cases, one hrst computes LHS and RHS as rational functions of kinematic data 
(in particular, they are functions of r), and then evaluate them on the same kinematics of 
(n+2) momenta that add up to zero. The statement of the double soft emission theorems is 
that the two are equal at the corresponding orders in r. It is this form of the soft theorems 
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that we check explicitly and prove in this section. 


A. Explicit Checks for the New Theorems 


For EMS, DBI and sGal theories, we write down the explicit form of the soft factors: 


= 7 E 


{ka-{p-q)Y 

^ti\ka-{p + q) + Tp-q 


+ ka-{p + q) + Tp-q 


2 n 




2 i^iK-{p + q) + rp-q 


t-3 n 

s(2) = I- y 


2 iAK-(p*q) + Tp-q 



■JT'f*(\-'ipq)(p-qf^ ■ 


(19) 

( 20 ) 
( 21 ) 


A fact that was not mentioned in Section I bnt crucial in the dehnition of and is 
that the operator Ja does not annihilate the prefactor (ka ■ {p+q) + Tp-q)~^ and so the above 
expressions literally mean that the prefactor comes in front of Ja- 

We computed pure scalar amplitudes in DBI theory analytically up to ten points, and 
checked the theorem at 0(r),O( t^),O( r^) for n = 4,6, and at for n = 8. We 

have also computed scalar amplitudes in EMS and sGal theory explicitly up to eight points, 
and conhrmed the complete soft scalar theorems for n = 4,6. 

In addition, we performed two important checks in DBI theory for amplitudes with exter¬ 
nal photons, (i) For the subleading soft scalar factors, the angular momentum of a photon 
has an additional spin part, Jq, photon = scalar+ -^a,spin, which acts on the polarization vectors 
('^aspin^a)^ “ ~ checked the leading and sub-leading theorems for 

the four-scalar-two-photon DBI amplitude, where in Jphoton is used for acting on the 
photons, (ii) We also checked the soft photon theorem for six-photon amplitudes in DBI. 
For partial amplitudes in YMS and NLSM theories, the soft operators read 



= 


kn • (p-q) + Tp-q ^ ki • {q-p) + rg 
kn ■ (p+q) + rp ■ q ki ■ (q+p) + rq ■ p j 
1 1 


T 


kn ■ {p+q) + rp-q 


PfiquJr + 


ki■(q+p) + rq- p 


qipPuJr 


( 22 ) 

(23) 


We computed NLSM amplitudes up to ten points, and checked the theorem at O(t^) 

for n = 4,6,8. We also conhrmed the theorem for six-scalar amplitudes in YMS theory. 
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B. CHY Representation and Soft Limits 


Before proceeding to the proofs, let us first review the CHY representation for amplitudes 
in these theories [18]. It is given by an integral over the moduli space of n-punctured Riemann 
spheres, with the locations specihed by holomorphic variables cr’s: 


/ n n p 

WdcraU'^(fa)In{W,k,...})=: / dfinin, (24) 

a=l a=l 

where the delta functions of fa '■= impose the so-called scattering equations [18], 

and the integrand In is some rational function that depends on the theory under consid¬ 
eration. The primes in the products denote redundancies in both the variables and the 
equations: for each product, one has to exclude three labels and compensate by a factor, 
e.g., n'a := (JidCTj^kCTk, (here aa,h '■= era - erf), and the result is independent of the 

choice. 

For scalar amplitudes of the theories we study in this paper, the integrand is a combi¬ 
nation of three basic building blocks [19]: (i) the Parke-Taylor factor for partial amplitudes 
with a given ordering, say with the canonical ordering (1,2,... ,n) 


{eri- (T2)(o-2 - 0-3) ■■■ {an 
(ii) the Pfaffian PfY„, and (iii) the reduced Pfaffian Vi'An, 
symmetric matrices are dehned by specifying their entries as 


-aiY 

where the two n x n anti- 



{Xn)ab := (1 - Sab) , {An)ab := (1 " Sab) 


an 


erb 


an. 


erb 


(26) 


Since the matrix An has co-rank 2 on the support of the scattering equations /a = 0, we 
dehne the invariant quantity Vi'An '■= ’j, where the minor is obtained by 

deleting the and rows and columns. With these building blocks, the integrands for 
scalar amplitudes of the theories we study are constructed as follows 


lf^^:={Vi'An)\ lY^^-.= ViXn{Vi'An)Y ■-= {ViXn)HVi'An)\ (27) 

2, ...,n):=C{l,2,...,n) (PfX)' , (1,2 ,... ,n) := C{1,2,... ,n) PfY„ Vi'An 


The integrals in (24) are localized on the a solutions to the scattering equations. With 
this formula the leading order of the above amplitudes under single and double soft emissions 
can be easily extracted by studying the behavior of the solutions for ex’s hrst. 
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For single soft emission, all the solutions are non-degenerate in the sense that {(Ja-CFb) ~ 
for any a, b. For double soft emission, however, apart from such non-degenerate solutions 
there exists a unique degenerate solution in which {cTn+i-<^n+ 2 ) ~ Given these we obtain the 
leading scaling in r of and of the building blocks on the solutions, which are summarized in 
Table I (we use “d” to denote the degenerate solution and “n” the non-degenerate ones). As a 


TABLE I. Leading Scaling of the Building Blocks 



single soft 

double soft (n) 

double soft (d) 

dfi 

T-l 

r-2 

T-l 

C 

rO 

rO 

T-1 

PfA 

rO 

rO 

r-1 

Pf'A 





consequence, the leading scaling of the formulas are summarized in Table If. Now it becomes 
obvious that double-soft behavior of these theories are special in that the contributions from 
non-degenerate solutions are suppressed: only the degenerate solution contributes to the 
hrst three and two orders for the two classes of theories respectively, which are exactly the 
orders we considered in the soft theorems. 


TABLE If. Leading Scaling of the Eormulas 



single soft 

double soft (n) 

double soft (d) 

sGal 




DBI 




EMS 



r-i 

NLSM 




YMS 

rO 

rO 

r-2 


C. Proof of the Leading-Order Double Soft Theorems 

Here we present the proof for the double soft scalar theorems at the leading order (see [20] 
for the proof at sub-leading order, and that for the double photon theorem.). As discussed 
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before it suffices to consider the degenerate solution only, for which we make the following 
transformation (we choose not to delete dan+idan+ 2 S(fn+i)S(fn+ 2 ))- 

^ 

O'n+l “ d ~ 2 ’ ~ 2 ' 

dan+ldan+2d{fn+l)d{fn+2) = -2 dpd^ 5{fn+l + fn+ 2 ) d{fn+l - fn+ 2 ) ■ (28) 


The basic idea here is to localize the ^-integral by 6{fn+i - fn+ 2 ), and regard the p-integral 
as a contour integration whose contour wraps the zeros of the equation fn+i + fn +2 = 0. Since 
^ = crn +2 - O'n+l ^ T, we Can further expand it as ^ Starting from the formula 

for Mn +2 this leads us to 


M, 


/ / 2m 


n+2 


ei 


yin ka-[p+q) q- n . n 
Zja=l p-aa 


In+2 + (sub-leading), 


(29) 


where the p-contour is now specihed by T,a=i ^ 0; evaluated on its unique 

solution ^ ELi 

To derive the leading-order soft theorems from (29), we need to expand J„ with respect 
to r. We hrst look at each building block, which behaves as 


C'(l,2,...,n + 2) = C'(l,2,...,n) 


On Oi 


(On - p) (-r ^ 1 ) (p - Oi) 




PfX 


n+2 ~ 


Pf'A 


n+2 


rp ■ q 


PfX„ + C>(T°), 


Pf'A„ + C>(r2). 


Combining these we obtain for scalar amplitudes in sGal, DBI and EMS (with m = 

^ (^2 / 1 \ 1-m / ^ ^ ^ \ m+S 


Mn+2 - 


fdpnln^ 2^^ 

dp. 




dn 


yn ka-(p+q) q- p . n 
Zja=l p^aa 

<f±. 

J 271'i 


(7k) (lr) 


.2m+2 


) 


/ o X , ka-(p-q) Y 

T{T^p-q)^ f dp j 


yn ka-(p+q) 
l-ia=l p-cTa 


+ 0{+^+^). 


(30) 

(31) 

(32) 


1 , 0 ,- 1 ) 


(33) 


Now we perform the p-integral by deforming the contour and use a residue theorem. Al¬ 
though there appears to be a simple pole at p = 00 , it is eliminated by an additional zero 
in the numerator due to momentum conservation. Thus we only encounter simple poles at 
p = Oa (a = 1,..., n), and the hnal result is 


Mn+2 = {r^p-qr 


(r ^ {ka-{p-q)f 
\4^i ka-{p + q) 


Mn + . 


(34) 
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Note that the soft operator in the bracket agrees with in (19) at 0{t), because the 
additional piece in is of higher order by momentum conservation, Ea=i ka ■ {p + q) = 
-2T{p.qy. This concludes our proof for the leading order soft theorem in these theories. 
Similarly, for scalar partial amplitudes in NLSM and YMS (with m = 0,-1) we obtain 


M(l,2,...,n + 2) 


f df^nln ^ 27 n 


dp iHp-q 






n dn 


p-aa 

(r^p-g)"* f dp 
2 J 27ri 






6 / 


, v^n kaip-q) 
dp L,a=l 


p-Ca 




Ea-1 {(^n-p) {p-f^l) 

^a-1 p-Ga 


+ C)(r"™+i), 


(35) 


Obviously there is no pole at p = oo, and we pick up two poles at p = cxi and p = cr„, which 
lead to the correct leading order soft theorem in these theories. 


M(l,2,...,n + 2) 


(^2p.g)m / ■ (p - q) ^ ki -jq-p) 
2 \kn-{p + q) ki-{q + p) 


M(l,2,...,n) + C>(r2™+i). (36) 


III. DISCUSSION 

Although not mentioned in Section I, the original motivation to search for subleading 
soft theorems in graviton emission in [1] was the connection among soft theorems and Ward 
identities of BMS symmetries [21, 22] at null inhnity [23]. Very recently, a Kac-Moody 
structure was found for four dimensional Yang-Mills also at null inhnity and two consecutive 
soft limits play an important role [24]. It is very tempting to suggest that the double soft 
theorems for scalars and the ones for photons could also have an interpretation as hidden 
symmetries. It is well known that in theories where single soft scalar limits vanish, double 
soft scalar factors carry information about non-linearly realized symmetries. It would be 
interesting to carefully explore the meaning of each of the soft theorems we found. 

As explained in Section IIB, the CHY integrand evaluated on different solutions to the 
scattering equations behaves in different ways, as summarized in Table II. There was a single 
degenerate solution (d) while all the others were non-degenerate (n). Denoting the leading 
double soft behavior of the integrand on a given solution as 0(r“d) or 0(r““), one hnds 
that the presence of soft factors is related to the difference 

A = On - CKd. (37) 
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We have seen that in all theories where A = 3 one hnds the presence of and 

5 ( 2 ) while in theories with A = 2 only are present. One can also compute A for 

the emission of other kind of particles, e.g., photons in DBI and EMS theory. In (15) we 
provided a very compact formula for S'O) of double soft photon emission and one might 
wonder whether there exists sub-leading soft factors. In both theories A = 1 which suggests 
that there is no other universal soft factor. Finding a clear explanation of the relation 
between A and the presence of soft factors is clearly one of the directions of future research. 

A related question is the derivation of the new soft theorems using BCFW techniques as 
was done originally in [1]. A clear hrst problem is that BCFW techniques present technical 
challenges when applied to scalar theories if no extra symmetries (such as supersymmetry) 
are present. However, we believe that A might somehow control the behavior for inhnity 
momenta of some generalization of BCFW appropriate for these theories. Very recently 
a complete analysis of applications of BCFW-like techniques was carried out in [25] and 
perhaps their approach can shed some light on this issue. 

We proved some of the soft theorems using the CHY representation, which appears to 
be particularly powerful for studying soft behavior of amplitudes. An important point to 
mention is that some of the CHY formulas here are still conjectural. This is why the checks 
done on explicitly known amplitudes are relevant. The still conjectural CHY representations 
have passed many non-trivial tests such as factorization. Therefore our proofs using them 
are very strong evidence for the validity of the leading and sub-leading soft theorems. 

Finally, let us mention an interesting analogy between the single soft theorem of a graviton 
(gluon) and the double soft theorem of theories of the hrst (second) class. Consider the 
single graviton soft factors and replace the polarization tensor by the product of two 
polarization vectors, i.e., Then 

0 ( 0 ) _ ^ 0 ( 1 ) _ ^ i^fika){ei^kn+l^pJa ) ^( 2 ) _ ^ ^ i^pkn+l^pJa^)^ 

- 2^ 1 “^gravity “ Zj U . U ’ gravity ~ r, 2^ ,, ' 


^ h . ■ k ^ “gravity h . ■ h ’ gravity n i i 

2=1 ^n+1 a=l ^n+1 ^ a=l ^n+1 

Comparing this form with the double soft factors for theories in the hrst class (7) it is easy 
to see that by identifying 


and Cl, soft graviton ^ Cl + C 2 (39) 

one hnds a striking similarity. The same analogy exists for theories of the second class and 
the single soft gluon emission. It would be interesting to sharpen this connection. 
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